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Abstract. We derive the TAP equations for the fermionic Ising spin glass. It is found 
that, just as in the non-fermionic model, the conditions for stability and for validity 
of the free energy are equivalent. We determine the breakdown of the paramagnetic 
phase. Numeric solutions of the fermionic TAP-equations at T = allowed to localize 
a first order transition between the spin glass phase and the paramagnetic phase at 
/i « 0.8. We computed at zero temperature the filHng factor i/(/Lt) and the distribution 
of internal fields. The saddle-point equations resulting from the calculation of the 
number of solutions to the TAP-equations were found to be much more complicated 
as in the non-fermionic case. 



Submitted to: J. Phys.: Condens. Matter 



Fermionic TAP-equations 



2 



Introduction 



Thouless, Anderson and Palmer derived local mean-field equations for the random Ising 



model \TE\. We generalize these TAP-equations to the four-state fermionic spin glass 
including a chemical potential /i. h is the number operator, Jij are gaussian distributed 
random interactions with variance J'^/N: 

ij i i 

It is believed 0, |1^, that the complete set of solutions to the TAP-equations is 
equivalent to the fully replica symmetry broken solution in the quantum-field theory 
0) 0' & El' so only known for half filling. We were able to solve the fermionic TAP- 
equations numerically on the T = axis and to determine the first order transition 
at /ic ~ 0.8 where the paramagnetic- and spin glass free energies become equal. The 
dependence of the filling factor on the chemical potential and the distribution of internal 
fields are given below. 



1. Fermionic TAP-equations 

The free energy corresponding to the non-fermionic TAP-equations is: 



^ = - E ^T'm. J-^^"^^"^^ - 1 E 4(1 - - 

i 

In order to calculate the corresponding expression for the generalized fermionic 
model, we extended the linked cluster diagrammatic theory by Horwitz and Callen ||^. 
For the sake of simplicity this rigorous and lengthy derivation is replaced by shorter and 
more intuitive arguments. 

The terms involving logarithms in equation correspond to the entropy of an 
ensemble of spins in the non-fermionic two state model with relative occupations 
and Uii, using rrii = rii-^ — Uii- Now we replace these terms with the entropy of the 
ensemble in the extended fermionic four state model. The relative occupations are 
denoted by nj^, n^j, n^o and setting rrii = rii^- — nj|, = rii^ + rii^ + 2nj|| and 
Qi = + nil- Then we account for the non-trivial occupation number of the magnetic 
states by replacing 1 — mf in the Onsager reaction field by qi — mf. After a final Legendre 
transformation the fermionic free energy reads 



Fermionic TAP-equations 



3 



ij ij i 

^ (2) 



2 

^ (l-g,)(l + T) / (1 - g~ )(1 + T) 



2 'V 2 

The system is characterized by the following 2 coupled TAP-equations. 

smh{(3Hi) 



m. 



cosh.{f3Hi) + cosh(/5/i) exp(-/5Xj) , . 

. cosh{(3Hi) ' 

cosh(/3ifj) + cosh(/9/i) exp{—f3Xi) 

where i/^ = /if* + J^j Jij'^j - P'^i J2j JiMj - and = f Jfji^j - ^J)- 
A third equation for the local filling factors Ui = 1 + (1 — gj) tanh(/3yu) follows 
from u = —d^Q. In the replicated quantum-field theory the corresponding equation 
V = 1 + (1 — g) tanh(/5/i) turns out to be invariant under an arbitrary number of replica 
symmetry breaking steps. 

2. Convergence and Stability 



Inspired by Plefka's work on the non-fermionic system and adopting his replacement 
Jij — * aJij we can rederive the fermionic free energy by a Taylor expansion as follows 



+ 



:i - + T) , / (l-g,)(l + T) ^ 



2 'V 2 



^ ^ (l-g.)(l-T) ^ /(I - - T) 



- /i ^(1 + (1 - g,) T) - ^ /^r^m, + ^ In 2 + 0{a'). 

i i 

The correct free energy is given by ^(a = 1). The diagrammatic expansion showed 
already that, provided the series actually converges, the terms of cubic order or higher 
are suppressed in the thermodynamic limit. We may now determine the radius of 
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convergence of ^(a), which, by a standard theorem of complex analysis, is equivalent 
to the radius of convergence of 

daQ{oi) = -^iYl "^ii^i^i)" = ~^Y1 J^i^i^i ~J6^ JijXij{a)- 

ij ij ij 

There, the susceptibility matrix x is defined as 

Xij{oi) = P{{cri(Tj)a - rriimj) = dh,dh^J^{h,x)- 

The radius of convergence is given by min(|a|), where the minimization is done 
over the values of a with eigenvalues in the inverse susceptibility matrix x~^- the 
non-fermionic case the equation {x~^)ij = {[dhdhJ-'{h)]~^)ij = 9m.(9m^.^(m) signifies that, 
when taking into account the special properties of the spectra of these random matrices 

the local stability of a TAP-solution implies the validity of the free energy at this 
point. In the fermionic model we have to deal with two different matrices to describe 
either the convergence of the linked cluster expansion {x~^)ij or the local stability of a 
given TAP-solution by 



A theorem by Pastur |]10|, ITl| used heavily by Plefka []T2| and resolvent- calculus 



can be applied to both cases to determine the limits of the support of the spectra. It 
is very interesting to note, that at the end of quite lengthy calculations both matrices 
lead exactly to the same set of conditions generalizing Plefka's convergence and stability 
conditions by 



{{q-m'f) <T' (5) 

i(g(l-g)) + 2(m2-m^) <T2. (6) 

The known non-fermionic limit is obtained by setting (ji = 1. We have thus linked 
the local stability of TAP-solutions to the finite support of the spectrum of random 
matrices in the thermodynamic limit. This gives a hint why the numerical search for 
these solutions is so difficult ^ for finite N the support of these spectra becomes 
unbounded, see e.g. for the exponential corrections to the semi-circle-law at finite 
N. This means, that for every solution the probability of having negative eigenvalues in 
the stability matrix is finite. But already one negative eigenvalue prevents finding this 
solution via minimization algorithms. Most solutions become thus unstable. 



3. Breakdown of the homogeneous paramagnetic solution 

The TAP-equations are readily solved numerically in the homogeneous paramagnetic 
phase, where they reduce to one single equation for q = (ji given for all i by 

1 + cosh(/3/i) exp(— ^^{^ q) ^ ^ 
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This is exactly the equation found in for q with a rephcated quantum-field 

theory approach. 

The second-order transitions between the paramagnetic and the spin glass phase 
are given by the intersection of the paramagnetic solutions with the stability conditions 
(equations (^) and (^). In the next section we find a first order transition on the T = 
axis at /Xc ~ 0.8. We expect a line of first order transitions linking this point with the 
tricritical point. The dotted line in the phase diagram (figure |l|) gives only a schematic 
behaviour of this line as the exact path is yet unknown. 
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Figure 1. Phase diagram obtained from tlie fermionic TAP-equations. The dotted line 
is a linear approximation th the first order transition line connecting the (calculated) 
tricritical point and the T = critical point. 



4. TAP-equations at T = 

At T = the TAP-equations reduce to 

rrii = Q{hi), qi = ml 

where hi = ^^^Jijini and VL{x) = 9(x — /i) — 0(/i — x) denotes a modified sign 
function. The energy corresponding to these solutions is simply 

fsG = - ^ Jijmirrij - /x ^ 

which has to be compared with the free energy of the paramagnetic solution 
fpM = — 2/i to find the first order transition. We were able to calculate numerically 
a huge number of spin glass solutions. We first note the interesting dependence of the 
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filling factor u on the chemical potential (see figure |^) . Unlike the discontinuous replica- 
symmetric and finite-step RSB solutions ^ the filling factor varies continously with 
/i in the vincinity of /i = 0. The numerical data for the increase of near fi = are 
compatible with power law fits oc x > 1, including exponential behavior. 
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Figure 2. Filling factor v(^[i) as a function of the chemical potential for systems of 
different size. 



From the dependence of the energy difference between the spin glass solutions 
and the paramagnetic solution (figure |) we can deduce a first order transition at 
= /^ci ~ 0.8. This critical value of can be viewed as a T = analogue of the 
Tci for the thermal first order transitions. 

Another very interesting feature shows up in the behaviour of the distribution of 
the local fields h^. When applying a chemical potential this probability density function 
is substantially modified. The "softgap" (see [T^) at /ij = splits up into two soft gaps 



at hi = — /i and hi = /i. Within the interval [— /i, /i] another peak emerges (see figures 
land I). 



5. Number of solutions 



The number of solutions of the TAP-equations (related to the so-called complexity) can 
be calculated adopting the procedure of Bray and Moore |I[ for finite T or following 
Roberts for T = 0. We were able to obtain in both cases fermionic generalizations of 
the saddle-point equations, which become extremely complicated due to the additional 
non-magnetic degrees of freedom. For example in the finite temperature case the 
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equations for the parameters g, q t], g, A and B read 




1 + B{k-{k\Y)' 



1 



= -(KA(tanh-VA) - kAA))/ - A - .fiq - q) 
J2 



(8) 



V 



-(A + B)- —{qVT^k^E)i 



A + B + 



2q 



J^q 



{{tanh-\\) - K\Af)j 



J2 



which can be viewed as an extension of the original equations given by Bray and 
Moore, but the average (•)/ means averaging by use of the following kernel: 



dn J d\g{x, K, \)S{f{x, K, X)) 



K 



1 (kA)2) 



(!-«:) 



X exp (t] {kXY + g k + ^ tanh ^(A) ^ (tanh ^(A) — kAA)' 



X 5 In 



'l-K 



A2 



J\q-q) 



X = itl/k is the reduced magnetisation, S an abbrevation for a very lengthy 
expression. These equations remain currently unsolved even numerically. The equations 
for T = are equally hard to treat. 



6. Outlook 



The numerical work presented here should now be accompanied by exact analytic RSB 
calculations for arbitrary /i. In order to reproduce the behavior of z/(/i) and of P{hi,fi) 
the generalizations of Robert's saddle-point equations |13| should be solved. It would 
be desirable to have Parisi's solution (infinite step RSB) for fi 0. The answer to 
one of the open questions might reveal the exact path of the first order transition line 
for T 7^ 0. If this turns out to be impossible, one should find more refined numeric 
algorithms, which allow to solve the fermionic TAP-equations for T ^ 0. 
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(a) iV = 75 




(b) N = 150 



Figure 3. Energy gap between spin glass solutions of the zero temperature TAP- 
equations and the paramagnetic solution for systems of different size. 
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Figure 4. Internal field distribution P{\hi\) for N = 150 and different values of fj,. 
The distribution P{hi) is symmetric. We used about 10000 points for each histogram. 




Figure 5. Scaling of P{hi) &t N = 50, N = 100, and N = 150 for = 0.6. 



